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We analyze the scattering properties of a periodic one-dimensional system at criticality represented 
by the so-called power-law banded random matrix model at the metal insulator transition. We focus 
on the scaling of Wigner delay times r and resonance widths V. We found that the typical values 
of r and T (calculated as the geometric mean) scale with the system size L as r typ cc L Dl and 



-(2-r> 2 ) 



where Di is the information dimension and D2 is the correlation dimension of 



eigenfunctions of the corresponding closed system. 



PACS numbers: 03.65.Nk, 71.30.+h, 72.15.Rn, 73.23.-b, 



I. INTRODUCTION 

The investigation of the disorder induced metal 
insulator transition (MIT) is an exciting and yet unsolved 
problem in condensed matter physics. During the past al- 
most fifty years of research many important features have 
been cleared out using analytical, numerical, and exper- 
imental techniques. Yet this problem continues to bring 
interesting and unexpected novelties to daylight i2*2i 

Many characteristics of disordered metals and those 
of insulators have been understood over the past. The 
spectral fluctuations of a disordered metal are well de- 
scribed by random matrix theory, while the fluctuations 
of a system in the insulating regime follow the Poisson 
statistics!^ The nature and the details of the MIT, on 
the other hand, still belong to the most intensively stud- 
ied problems. We know already that the spectral fluctu- 
ations of a system right at the transition have properties 
merged from the two extremes^ However, relations be- 
tween the spectral and eigenfunction fluctuations show 
that there is an intimate connection between the twc& 
which is related to anomalous wave packet spreading and 
diffusion. ' 

Recently, several works have been devoted to deepen 
our understanding of the scattering properties of disor- 
dered systems by analyzing the distribution of resonance 
widths and Wigner delay times. 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 : 16 ! 17 Both 
distribution functions have been shown to be closely re- 
lated to the properties of the corresponding closed sys- 
tem, i.e., the fractality of the eigenstates and the critical 
features of the MIT. In this respect, detailed analysis 
have been performed for the three-dimensional (3D) An- 
derson modeU^ and also for the power law band random 
matrix (PBRM) model. 1 ?! 1 ?! 21 ?! 21 ' 22 ! 23 The PBRM model, 
which describes a one-dimensional (ID) sample with ran- 
dom long-range hopping, has been found to provide many 
of the features of the localization-delocalization transi- 
tion present in the 3D Anderson model. 

In the present work we look at the PBRM model more 



close ly to see how general are the findings presented in 
Refs.im and[l]l Here, we find good scaling of the dis- 
tribution functions of resonance widths and Wigner delay 
times. Moreover, we state new scaling relations of their 
typical values (i.e., their geometric mean) determined by 
the information dimension, D\, and the correlation di- 
mension, Z?2 , of the eigenfunctions of the corresponding 
closed system. We note that D2 also governs the spread- 
ing of wave packets in systems at the MITii 

Finally, we have to mention that in recent microwave 
experiments the study of systems in which the resonance 
states resemble the characteristics of those of our interest 
was reported^ Therefore the direct verification of our 
results is expected to be available in the near future. 

The organization of this paper is as follows. In the 
next section we describe the model we use and define the 
scattering setup. Section III is devoted to the analysis 
of the coupling between sample and lead. In Sec. IV 
we present our results for the distribution functions of 
resonance widths and Wigner delay times. Finally, Sec. 
V is left for conclusions. 



II. MODEL 

The isolated sample of length L is represented by an 
L x L real symmetric matrix whose entries are ran- 
domly drawn from a normal distribution with zero mean, 
(Hij) = 0, and a variance depending on the distance of 
the matrix element from the diagonal as 
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2l+[sin(7r|i-i|/L)/(7r&/L)] 
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where b and a are parameters. In order to reduce finite 
size effects this expression already incorporates periodic 
boundary conditions, where Eq. Q is known as the peri- 
odic PBRM model. Nevertheless, for sites far away from 
each other and the border, i.e., 1 <§; \i — j\ <C N, the 
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variance decays with a power law 
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Field-theoretical considcrationsiS'SLSi and detailed nu- 
merical investigations^*^ verified that the PBRM model 
undergoes a transition at a = 1 from localized states for 
a > 1 to dclocalized states for a < 1. This transition 
shows all the key features of the Anderson MIT, includ- 
ing multifractality of eigcnfunctions and non-trivial spec- 
tral statistics at the critical point. At the center of the 
spectral-band a theoretical estimation for the multifrac- 
tal dimensions, D q , of the eigenfunction ip(r) gives^ 

n = Ubf(q-l/2)[^(q-l)T(q-l)]-\ 6«1 
q \ l-q(2irb)-\ 6»1 

where T is the Gamma function. D q is defined through 
the Renyi entropy^ of the participation numbe»22i2i2Si2i 



1Z„ = In J\f Q oc D„ In L 

1 — q 



(4) 



where the participation number is defined as 

K=(J Wr)| 29 dr) (5) 

showing a nontrivial scaling with respect to the linear size 
of the system L. Thus model possesses a line of crit- 
ical points b G (0, oo), where the multifractal dimensions 
D q change with b. The value of parameter b describes the 
strength of disorder. In other words it is related to the 
classical conductance through the relation^ g c = 4/36, 
where (3 stands for the global symmetry of the system: 
(3 = 1 describes the presence and (3 = 2 the absence of 
time reversal symmetry. 

Among all dimensions, the information dimension D\ 
and the correlation dimension Di play a prominent 
roleiSS However, D\ is obtained as the q — ► 1 limit of 
Eq. being the scaling of the Shannon entropy^ 

We have to mention that recently a phenomenological 
formula was obtained in Ref. for the exponent D2 as 
a function of b (and in this way as a function of g c ): 
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where a is a fitting parameter. For the present model we 
found that a m 2.85 reproduces well the numerical D2 
extracted from Nir^ see Eq. ©. Similarly Eq. ® fits 
well the exponent D x with a « 5.6. See Fig. ^ 

We turn the isolated system to a scattering one by at- 
taching one semi-infinite single channel lead to it. The 
lead is described by the following ID tight-binding Hamil- 
tonian: 



H 



lead 
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FIG. 1: D\ and D2 as a function of b. Lines correspond to 
Eq. © with a w 5.6 (dashed) and a m 2.85 (dot-dashed). 



Using standard methods^ one can write the scattering 
matrix (a single complex number for one attached lead) 
in the formic 



S(E) 



JHE) 



l-2iirW T {El-Hctf^W, (8) 



where 1 is an L x L unit matrix and 7i c ff is an effective 
non-hcrmitian Hamiltonian given by 



n eS = H-inWW- 



(9) 



Here, W is an L x 1 vector with elements W n = woS nno , 
where wq is the coupling strength between sample and 
lead and no is the site at which the lead is attached. All 
our calculations take place in an energy window close to 
the band center [E = 0). Then, the Wigner delay time 
is given by 14 i 16 i 31 



t(E = 0) 



dE 



-2ImTr(S - H eS ) 



E=0 



\E=0 



(10) 



In addition to the delay times, which captures the 
time-dependent aspects of quantum scattering, the poles 
of the scattering matrix arc also of great relevance^ 
They determine the conductance fluctuations of a quan- 
tum dot in the Coulomb blockade regime^ or the cur- 
rent relaxation^ The poles of the scattering matrix show 
up as resonances, which are the complex eigenvalues 
£ n = E n — «r n /2 of 7i c ff. E n and T n are the position 
and width of the nth resonance, respectively. Of course 
for wq — the real part of the poles, E n , correspond 
to the eigenvalues of the closed system. The r„s, on 
the other hand, are related to the lifetime of the reso- 
nances through r n = l/r„. Thus, we may expect close 
correspondence between Wigner delay times and the res- 
onance widths which, in fact, we calculate independently 
from each other. Bellow we use matrices of sizes varying 
from L = 50 up to L = 800 to compute probability dis- 
tribution functions of scattering phases Wigner delay 
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FIG. 2: corresponding to wo where \{S)\ takes its min- 

imum value (perfect coupling condition) for (a) b = 0.01 and 
0.04, (b) b = 0.1, (c) b = 0.4, (d) b = 1, (e) b = 4, and (f) 
b = 10. We use (a) too = 0.6, (b) w = 0.65, (c) w = 0.9, (d) 
too = 1.55, and (e) wo — 3.4. In (f) too = 5.1, wo = 5.43, and 
too = 5.5 for L — 50, L — 100, and L > 200, respectively. 



times t, and resonance widths T. For statistical process- 
ing a large number of disorder realizations is used. In the 
case of V{&) and V(t) we have at least 100000 data val- 
ues. For the statistics of T we use the eigenvalues around 
E n ~ 0, about one tenth of the total spectra. Here, at 
least 250000 data values are used. 
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FIG. 3: (color online) lnP(lnr/r typ ) as a function of lnr/r typ 
for (a) b = 0.1, (b) 6 = 1, and (c) b = 10. In the insets the 
scaling r typ cx L~ M is shown together with a linear fitting. 
r typ = exp(lnr). We found (a) fi = -0.318 ± 0.027, (b) 
H = -0.874 ± 0.006, and (c) (j, = -0.997 ± 0.016. 



III. PERFECT COUPLING 

In this section we define more precisely the way the 
lead is attached to the system under the condition of 
perfect coupling. We should first note that the lead 
is attached to any site (chosen at random) within the 
sample and that the ensemble average also implies av- 
erage over the position of the lead. The value of the 
coupling strength, wq, is defined by analyzing the distri- 
bution of the scattering phases V{<&). See also Ref. 0- 
In Fig. [3 we present V($>) for various values of b and 
L. We use (here and bellow) the coupling strengths 
wq such that |(S(i/;o))| ~ 0, where the corresponding 
setup is associated with perfect coupling between sam- 
ple and lead. (S) is the ensemble average of the S 
matrix. Notice that for b > 1 we recover the uni- 



form distribution. Moreover, in this case, the expres- 
sion£'P(<I>) = l/27r(7+^7 2 - lcos($)) holds (not shown 
here); where 7 = (1 + |(5)| 2 )/(1-|(5)| 2 ) and < (S) < 1. 



IV. DISTRIBUTION OF DELAY TIMES AND 
RESONANCE WIDTHS 

We first have to mention that moments of the the in- 
verse Wigner delay time have been (numerically and an- 
alytically) shown to be related to the multifractal dimen- 
sions of cigenfunctions as&2iiS 

(t _ «) =L- qD "+\ (11) 
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FIG. 4: (color online) In V int (In r/r typ ) as a function of 
In r/r typ for (a) b = 0.1, (b) 6=1, and (c) b = 10. The 
thin dashed line in (c) with slope -3/2 is plotted to guide the 
eye. 



although the validity of this relation depends strongly 
whether the lead is attached to a site inside or at the 
border of the sample. However, in the case of the peri- 
odic PBRM model, studied here, all sites are generic and 
Eq. I|ll|) holds well (not shown here). 

Now, let us turn to the main results of our present in- 
vestigation. We have looked at the distribution of Wigner 
delay times as well as that of the resonance widths a lit- 
tle bit differently, compared to previous works. First we 
tried to obtain a universal (i.e., L independent) form of 
the distribution function 'P(lnTA), where A oc 1/ (r) ~ 
1/L. However, for small b the relation (r) ~ L does not 
hold and we have found a breakdown of the scaling of 
■p(hiTA) with L, in contrast to the expectation accord- 
ing to Ref. El 

Instead of (r) we decided to look at the typical value 
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FIG. 5: (i as a function of b. n is extracted from the scaling 
r typ tx L" M , see Fig. 01 The dashed line is -Di, where Di 
was obtained from Eq. JfjJ with a « 5.6. 



of r: r typ = exp (lnr). We found that T>{h\ r/r typ ) shows 
a good scaling with L and a tendency to a universal dis- 
tribution (for each b) as L — > oo, sec Fig. |3 In Fig. 0] we 
also plot V iat (In r/r typ ) to see the behavior of the tails of 
V (In r/r typ ), where 



V{x')dx' . 



(12) 



Here, only for b = 10 we see a clear decay of the form 
V IDt (T/T typ ) ~ (r/r typ )- 3 / 2 . Notice that the value of the 
exponent is the same as the one found in Ref. ^] for 
PmtOV (r» ~ (t/ (t))- 3 / 2 in the case of the 3D Ander- 
son model. 

In view of Eq. |JTJ we may estimate the scaling 
behavior of r typ . Using a simple relation^ (lnr) = 
lim„ 



i9 _ (<T 9 ) - l)/q, we obta 



oc L 



(13) 



This is contrasted with the numerical results presented 
in Fig. \5\ where we plot /i as a function of b. Here, fj, is 
extracted from the scaling r t; 



x 



(see the insets in 
Fig. OJl. The data agree well with Eq. ((131) . 

In the case of resonance widths, we observe an excel- 
lent scaling of V(lnT/r typ ) for all values of b and the ap- 
proach to a universal distribution (for each b) as L — > oo, 
see Fig. El Here, T typ = exp (In T) . Also, from Fig. □ we 
observe that as L — ► oo, the integrated distribution func- 
tion follows a power law decay V iDt (T /T typ ) ~ (r/r typ ) _1 
for all values of b in complete agreement with Refs. ITil 
andll7i 

Finally, in Fig. we plot vasa function of b where v is 
extracted from the scaling r typ oc L~ v (see the insets in 
Fig. . We found that v is well approximated by 2 — D 2 - 

Notice that in the limit b — > oo, where v — ► 1, we 
recover the full random matrix limit in which case we also 
have r typ oc L oc 1/A and r typ oc oc A. This means, 
for instance, that as b -> oo, P(r/r typ ) -> V(TA), where 
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FIG. 6: (color online) ln7 3 (lnr/r typ ) as a function of 
lnr/r ,yp for (a) b = 0.1, (b) 6 = 1, and (c) b = 10. In the 
insets the scaling F typ oc L~ v is shown together with a linear 
fitting. r typ = exp(lnr). We found (a) v = 1.728 ± 0.005, 
(b) v = 1.219 ± 0.007, and (c) v = 1.048 ± 0.008. 



the later scaling was found in Ref.^|for the 3D Anderson 
model. 

In Figs. [S] and El we see that the two quantities, r typ 
and r typ , behave very similarly and in the limit of weak 
multifractality, b ^> 1, we nicely recover the features of 
the 3D Anderson transition^ Moreover, this similarity 
calls for a simple scaling law 
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FIG. 7: (color online) ln-p int (lnr/r typ ) as a function of 
lnr/r ,yp for (a) b = 0.1, (b) 6 = 1, and (c) b = 10. The 
thin dashed lines with slope -1 are plotted to guide the eye. 
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(14) 




V. CONCLUSIONS 

We have presented an analysis of the scattering prop- 
erties of the PBRM model with one lead attached to it. 



FIG. 8: v as a function of b. v is extracted from the scaling 
F typ oc L~" , see Fig.© The dashed line is 2 - D 2 , where D 2 
was obtained from Eq. JHJ with a ~ 2.85. 
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Especially, the system with large conductance (weak mul- 
tifractality) shows remarkable correspondence with the 
properties of the standard 3D Anderson model at criti- 
cality. This is in accordance with earlier findings. We 
have also found a novel relation between the scaling of 
the Shannon entropy, described by Di, and the typical 
values of the Wigner delay times; as well as a relation 
between the scaling exponent of the participation num- 
ber A/2, described by D 2 , and the resonance widths. On 
the other hand the decay of their integrated distribution 
functions behave just as they do for the Anderson tran- 
sition in three dimensions. 
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